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We report on the computation of a class of integrals that appear when integrating the so-called iterated singly- 
unresolved approximate cross section of the NNLO subtraction scheme of Refs. [US], over the factorised phase 
space of unresolved partons. The integrated approximate cross section itself can be written as the product of 
an insertion operator (in colour space) times the Born cross section. We give selected results for the insertion 
operator for processes with two and three hard partons in the final state. 
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1. Introduction 

High energy particle collisions frequently lead 
to final states with hadronic jets. Jet observables 
can be ideal for precision studies, since their large 
production cross sections allow them to be mea- 
sured with high statistical accuracy. Examples in- 
clude: the determination of the strong coupling, 
as J from jet rates and event shapes in electron- 
positron annihilation; the measurement of gluon 
parton distribution functions (and also as) in 
deep inelastic lepton-hadron scattering into two 
plus one jets; the determination of parton distri- 
butions in hadron-hadron scattering from single 
jet inclusive production and vector boson plus jet 
production. The relevant observables are often 
measured with experimental precision of a few per 
cent or better, thus theoretical predictions with 
the same level of accuracy are necessary. This 
usually requires the computation of next-to-next- 
to-leading order (NNLO) corrections in perturba- 
tive QCD. 

One of the main bottlenecks of the straightfor- 
ward application of QCD perturbation theory at 
NNLO is the following. Due to the presence of 
infrared (IR) singularities, the finite higher-order 
corrections are sums of pieces which are sepa- 
rately divergent in d = 4 spacetime dimensions. 
These IR singularities must be regularised and 
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cancelled before any numerical computation may 
be performed. The most common approach for 
handling IR singularities at next-to-leading or- 
der (NLO) accuracy is the subtraction method. 
After regularising all contributions by dimen- 
sional regularisation in d = 4 — 2e dimensions, 
one builds subtraction terms that simultaneously 
cancel both the kinematical singularities in real- 
emission phase space integrals and the explicit 
e-poles in one-loop virtual corrections [5]. Since 
the IR singularity structure of QCD amplitudes is 
universal, such subtraction terms can be defined 
in a general, i.e. process and observable indepen- 
dent fashion. 

In recent years, severe efforts have been made 
to extend the subtraction method to NNLO ac- 
curacy [TMIMT^. and this has proved to be a 
challenging problem. In broad terms, any pro- 
posed subtraction scheme must address two quite 
distinct difficulties. First, one must define sub- 
traction terms that properly regularise the real- 
emission phase space integrals. In a rigorous 
mathematical sense, the cancellation of kinemat- 
ical singularities must be local, i.e. the subtrac- 
tion terms must have the same pointwise singular 
behaviour in d dimensions, as the real-emission 
pieces. Second, one must combine the integrated 
form of these counterterms with the virtual con- 
tributions, cancelling the IR divergences of the 
loop matrix elements. Again, from the point 
of view of mathematical rigour, this cancellation 



must be local, i.e. it must happen pointwise in 
phase space. Practically, the full locality of the 
subtraction scheme is also important to ensure 
good numerical efficiency of the algorithm. Fi- 
nally, the construction should be universal (pro- 
cess and observable independent), otherwise te- 
dious adaptation of the algorithm to every spe- 
cific problem becomes necessary. 

The construction of a general subtraction 
scheme with fully local counterterms requires a 
lot of careful analytical calculations. In order to 
spare some work, various approaches have been 
proposed to NNLO calculations sacrificing either 
full locality or generality for easier analytical (but 
more cumbersome numerical) treatment. For ex- 
ample, in the antenna subtraction scheme [9^1 2j 
azimuthal correlations in gluon splitting are not 
reproduced (and the cancellation of e-poles in the 
real- virtual contributions is also nonlocal). As 
a consequence, actual numerical computation of 
total rates jl3f[T6] and event shapes [T7H2T| in 



e"'"e~ ^^ 3 jet production in the antenna scheme 
requires the use of an auxiliary phase space slic- 
ing. Another option is to develop dedicated 
subtraction schemes that are applicable only to 
some specific processes, such as the production 
of colourless final states in hadron-hadron colli- 
sions |22f[24] . Then, one may even propose to 
dispense with the subtraction method altogether, 
and adopt a strategy such as sector decomposi- 
tion (see e.g. [IS] and references therein). 

Nevertheless, it is possible to address all sub- 
tleties, and define completely local counterterms 
for real radiation, as done in full detail in Refs. 
[iHl] for the case of colourless initial states. 
(Work towards an extension to hadron-initiated 
processes is presented in Ref. [26].) Then, to 
make this subtraction scheme an effective tool, 
one must compute the integrals of all subtraction 
terms over the phase spaces of unresolved emis- 
sions. In this contribution, we briefly review the 
methods used to compute the integrated subtrac- 
tion terms, and report on the evaluation of the in- 
tegrated iterated singly-unresolved cross section, 
the term labelled A12 , which appears in Eq. ([5]) . 
However, we begin by recalling the essentials of 
the subtraction scheme of Refs. [IH1]- 



2. Subtraction at NNLO 

Consider the NNLO correction to a generic m- 
jet observable. 
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The three contributions on the right hand side 
are separately divergent in d = 4 dimensions, but 
their sum is finite for IR safe observables. To ob- 
tain the finite NNLO correction, we first continue 
analytically all integrals to d = 4 — 2e dimensions 
and then rewrite Eq. ([T]) as 
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that is a sum of three finite integrals where the 
integrands. 
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are integrable in four dimensions by construc- 
tion. The approximate cross sections da^^'g ^ 

RR A 

and da„j_|_'2 ^ regularise the doubly- and singly- 
unresolved limits of the real-emission piece, 
do'm+2 respectively. The double subtraction due 
to the overlap of these two terms is compensated 

RR A 

by da„_,l'2 " . These terms are given explicitly in 



Ref. [3]. Finally, dcr^+f^ and ( /^ da^+a"^' ) ^' 

regularise the singly-unresolved limits of dcrj^^^ 

and J^ 'io'm+2 ^ respectively. They are given ex- 
phcitly in Ref. [1]. 

The construction of each approximate cross 
section in Eqs. ©-(HI) is based on the known 
and universal IR limits of tree level and one-loop 
squared matrix elements, and proceeds in two 
steps. First, the IR factorisation formulae are 
written in such a way that their complicated over- 
lap structure can be disentangled ("matching of 
limits" ) jTll??] . Second, we define "extensions" of 
the formulae, so that they are unambiguously de- 
fined away from the strict IR limits |2H1]- These 
extensions are defined by the use of various mo- 
mentum mappings that map a set of rn -I- 1 or 
771 + 2 momenta into a set of m momenta. 
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such that (i) the delicate structure of cancella- 
tions among the matched limit formulae in vari- 
ous limits is respected (ii) exact momentum con- 
servation is implemented, and (iii) the original 
m + 1 or m + 2 particle phase space factorises 
exactly into the product of an m particle phase 
space and a one- or two-particle phase space mea- 
sure, 

d0m+i({p}m+i;Q) = d0„({p},„;Q)[dpi^„i] (7) 
and 

d0m+2({p},„+2; Q) = d0„({p},„; Q)[dp2,m] ■ (8) 

To finish the definition of the scheme, one must 
compute once and for all the one- and two-particle 
integrals, denoted formally as J^ and J^, appear- 
ing in Eqs. Q and ([5]). We discuss this next. 

3. Integrating the counterterms 

The actual computation of the integrated coun- 
terterms leads to a large number of multi- 
dimensional integrals. The ultimate goal is to 
find the analytical form of the coefficients of a 
Laurent-expansion (in e) of these integrals, which 
turns out to be a rather a tedious job. In or- 
der to compute these coefficients as efficiently as 
possible, we have explored several methods. 



First, it is possible to extend the method 
of integration-by-parts identities and solving of 
differential equations, developed for computing 
multi-loop Feynman integrals [IHlll^: to the rel- 
evant phase space integrations j30]. This method 
yields e-expansions with fully analytical coef- 
ficients, with the final results being expressed 
in terms of two-dimensional harmonic polyloga- 
rithms (after a suitable basis extension, see Ref. 
[3(J' for details). This approach was used suc- 
cessfully to compute a class of singly-unresolved 
integrals [30] . 

Second, the phase space integrals that arise can 
be computed via the method of Mellin-Barnes 
(MB) representations [5TH55] . Here we obtain 
the e-expansion coefficients in terms of complex 
contour integrals over F-functions. Performing 
these integrals by the use of the residue theo- 
rem, a representation in terms of harmonic sums 
is obtained. In many cases, the sums can be 
evaluated in a closed form, yielding an analyti- 
cal result. In some instances however, we find 
multi-dimensional MB integrals that are very dif- 
ficult to compute fully analytically. Nevertheless, 
in these situations a direct numerical evaluation 
of the appropriate MB representations provides a 
fast and reliable way to obtain final results with 
small numerical uncertainties. We stress that for 
phenomenological applications, this is all that is 
required, since the relative numerical uncertainty 
associated with phase space integrations is gen- 
erally much larger than that of the integrated 
counterterms. We have used the MB method to 
compute all singly-unresolved integrals [S^, and 
very recently to evaluate all two-particle integrals 
appearing in J^ dcr^_(_'2 ^^ as well [35] (see below). 

Finally, the method of iterated sector decom- 
position [3S] can also be used to calculate the 
integrals we encounter |:36 . Sector decomposi- 
tion produces a representation of the e-expansion 
where the coefficients are given in terms of 
(mostly quite cumbersome) finite integrals over 
the unit hypercube. The analytical evaluation of 
these integrals is not feasible except for the sim- 
plest cases, and we find that in most instances, 
the MB method provides an integral representa- 
tion for the expansion coefficients which is better 



suited for direct numerical integration. Neverthe- 
less, this method is simple to implement and we 
employed it to numerically cross check all our fi- 
nal results. 

4. Results 

In previous publications, all one-particle inte- 
grals, denoted formally by /^ in Eqs. dH) and dS]) 
have been evaluated with the methods just dis- 
cussed [SnUSHIMI- Here we report on the com- 
putation of the /„ d(T^ I 2 ^^ two-particle inte- 
grated counterterm '351 . To begin, we recall that 
the iterated singly-unresolved subtraction term, 
da^^2 ^^ can be written symbolically as [3] 
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where ^i2|A^m+2p further evaluates to a sum of 
several terms, as spelt out in detail in Ref. [3]. 
Each of these terms is essentially a product of 
kinematical factors times an m-parton factorised 
matrix element, the later being independent of 
the variables in the factorised phase space mea- 
sure [dpi_m][dpi,m+i]- Thus, it is possible to eval- 
uate the integrals of the kinematical factors over 
the factorised phase spaces, without reference to 
the specific process or observable. The final re- 
sult, after summation over unobserved flavours, 
can be written in the form of an insertion opera- 
tor (in colour space) times the Born cross section 
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Here the insertion operator has three terms ac- 
cording to the possible colour structures 
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with fi denoting flavours, Cq = Cp = T^, Cg = 
Ca = Tg, and S, = r^^- In Eq. ^, the de- 
pendence of the functions Cj^2 f i ^tc. on the kine- 
matics is suppressed for the sake of simplicity. 
These functions in turn are given as the following 
combinations of flavour summed integrated coun- 
terterms 
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kinematical dependence. We find 
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On the right hand sides of these equations, the 
flavour summed integrated counterterms them- 
selves are sums of various integrated subtraction 
terms. Their precise definitions are somewhat 
long, and will not be reported here. 

Instead, we present results for the insertion op- 
erator in Eq. (|11|) for processes with two and three 
partons in the final state. Consider first the pro- 
cess e+e^ -^ 2 jets. The corresponding Born ma- 
trix element is jA^j (I5, 2q)p, i.e. the quark car- 
ries the label 1 and antiquark label 2. The colour 
algebra is trivial fSl , and all properly scaled kine- 
matical invariants are equal to one, i.e. I12 {{p}2) 
does not depend on the kinematics. Introducing 
the notation [37j x = Ca/Cf, y — Tr/Cf, we 
find 
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where yo S (0, 1] and D'q > 2 and integer are fixed 
parameters that enter the precise definition of the 
subtraction terms. Finally, the function Y,{z,N) 
is defined as follows 
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fe=i 
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Next, consider e^e —5- 3 jet production. The 
Born matrix element is \M.]^ {lq,2q,3g)\'^, i.e. the 
quark carries label 1, the antiquark label 2 and 
the gluon carries label 3. The colour algebra is 
still trivial [S], but ^12 ({pja) now carries genuine 
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In this equation yik = 2pi -pk/Q^, with Q^ being 
the total incoming momentum. Higher order ex- 
pansion coefficients (in e) in Eq. (J19p are already 
quite cumbersome, and will be given elsewhere in 
the form of a computer program [35j . 

5. Conclusions 

In this contribution, we have reported 
the computation of the integrated iterated 
singly-unresolved approximate cross section, 
J^ dCT„^'2 " , that appears in the NNLO subtrac- 



tion scheme of Refs. [2li4l. The phase space 
integrals appearing in the computation can be 
evaluated once and for all, and their knowledge 
(as Laurent expansions in e, up to and including 
0(e") terms) is necessary to make the subtraction 
scheme an effective tool. All phase space inte- 
grals were evaluated with two separate methods. 
First, we used the method of Mellin-Barnes rep- 
resentations with harmonic summation to obtain 
analytical results where feasible. In some cases, 
obtaining complete analytical answers is very dif- 
ficult. In these situations, we integrated the MB 
representations of the expansion coefficients di- 
rectly, and this provides a fast and reliable way 
to obtain final results with small numerical un- 
certainties. Second, all integrals were computed 
with the method of sector decomposition as well, 
providing useful numerical checks. 

With the evaluation of J^ dcr,^^'2 ^^ , the last 
task in finishing the definition of the subtrac- 
tion scheme is the computation of the integrated 



doubly-unresolved approximate cross section, i.e. 
the term labelled A2 in Eq. ([S|). The analytical 
structure and complexity of the integrals that ap- 
pear in this final piece are essentially the same as 
the integrals considered here. Therefore, the tech- 
niques outlined above will also be applicable to 
the computation of this remaining contribution. 
This is work was supported in part by the 
Deutsche Forschungsgemeinschaft in SFB/TR 9, 
the Helmholtz Gemeinschaft under contract VH- 
NG-105, the Hungarian Scientific Research Fund 
grand OTKA K-60432. 

REFERENCES 

1. G. Somogyi, Z. Trocsanyi, and V. Del Duca, 
JHEP 06, 024 (2005) , arXiv:hep-ph/0502226, 

2. G. Somogyi and Z. Trocsanyi, (2006), 
'arXiv:hep-ph/0609041' 

3. G. Somogyi, Z. Trocsanyi, and V. Del Duca, 
JHEP 01, 070 (2007) , | arXwdi ep-ph/0609042 , 

4. G. Somogyi and Z. Trocsanyi, JHEP 01, 052 
(2007), arXiv:hep -ph/0609043, 

5. S. Catani and M. H. Seymour, Nucl. Phys. 
B485, 291 (1997),"arXiv:hep-ph/9605323' 

6. S. Weinzierl, JHEP 03, 062 (2003), 



arXiv:hep-ph/03 02180 

JHEP 



07, 052 (2003), 



7. S. Weinzierl, 
arXiv:hep-pli/0306248, 

8. S. Frixi one and M. Grazzini, JHEP 06, 010 
(2005), arXiv:hep-ph/0411399' 

9. A. Gehrmann-De Ridder, T. Gehrmann, and 
E. W. N. Glover, JHEP 09, 056 (2005), 



|arXiv:hep-ph/0505111 

A. Daleo, T. Gehrmann, 



and D. Maitre, 



10. 

JHEP 04, 016 (2007) , jarXiv:hep-ph/0612257 

11. A. Daleo, A. Gehrmann-De Ridder, 
T. Gehrmann, and G. Luisoni, JHEP 
01, 118 (^2010). [arXi^:0912.0374i 

12. E. W. N. Glover and J. Pires, (2010), 
larXiv: 100372 824 

13. A. Gehrmann-De Ridder, T. Gehrmann, 
E. W. N. Glover, and G. Heinrich, JHEP 11, 
058 (2007), arXiv:0710.0346, 

14. A. Gehrmann-De Ridder, T. Gehrmann, 
E. W. N. Glover, and G. Heinrich, Phys. Rev. 
Lett. 100, 172001 (2008), arXiv:0802.0813, 

15. S. Weinzierl, Phys. Rev. Lett. 101, 162001 



(2008) , ■arXiv:0807.3241' 

16. S. Weinzierl, JHEP 07, 009 (2009), 
iarXiv:0904. 11451 

17. A. Gehrmann-De Ridder, T. Gehrmann, 
E. W. N. Glover, and G. Heinrich, Phys. Rev. 
Lett. 99, 132002 l'2007'). [arXi^:0707. 12851 

18. A. Gehrmann-De Ridder, T. Gehrmann, 
E. W. N. G lover, and G. Hei nrich, JHEP 12, 
094 (20071 ■ larXiv:071 1.471 II 

19. A. Gehrmann-De Ridder, T. Gehrmann, 
E. W. N. Glover, and G. Heinrich, JHEP 05, 
106 (2009), ^^0903.4658, 

20. S. W einzierl, JHEP 06, 041 (2009), 
[ar^:0904.1077, 

21. S. Weinzierl, Phys. Rev. D80, 094018 (2009), 
arXiv:0909.5056, 

22. S. Catani and M. Grazzini, Phys. Rev. Lett. 



23. 



24. 



98, 222002 (2007), |arXiv:hep-ph/0703012 
S. Catani, L. Cieri, G. Ferrera, D. de Flo- 
rian, and M. Grazzini, Phys. Rev. Lett. 103, 
082001 (20091. iarXiv:0903.2120 , 
S. Catani, G. Ferrera, and M. Grazzini, 
JHEP 05, 006 (2010), arXiv:1002.3115, 

25. G. Heinrich, Int. J. Mod. Phys. A23, 1457 
(2008),'arXiv:0803.4177, 

26. G. Somogyi, JHEP 05, 016 (2009), 
larXiv:0903.12181 

27. Z. Nagy, G. Somogyi, and Z. Trocsanyi, 
(2007), arXiv:hep-ph/0702273, 

28. A. V. Kotikov, Phys. Lett. B267, 123 (1991). 

29. E. Remiddi, Nuovo Cim. AllO, 1435 (1997), 
arXiv:hep-th/ 9711188. 

30. U. Aglietti, V. Del Duca, C. Duhr, G. Somo- 
gyi, and Z . Trocs anyi, JHEP 09, 107 (2008), 
[arXiv:080735il 

31. V. A. Smirnov, Phys. Lett. B460, 397 (1999), 
arXiv:hep-ph/9905323, 

32. J. B. Tausk, Phys. Lett. B469, 225 (1999), 



arXiv:hep- ph/9909506, 

33. V. A. Smirnov, Springer Tracts Mod. Phys. 
211, 1 (2004). 

34. P. Bolzoni, S.-O. Moch, 
Z. Trocsanyi, JHEP 
[arSy:0905.4390, 

35. P. Bolzoni, G. Somogyi, and Z. Trocsanyi, 
(2010), in preparation. 

36. G. So mogyi and Z. Troc sanyi, JHEP 08, 042 
(20081 ■ rarXiv:0807.05"09l 



G. Somogyi, and 
08, 079 (2009), 



37. Z. Nagy and Z. Trocsanyi, Phys. Rev. D59, 
014020 (1999), arXiv:hep-ph/9806317, 



